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1. Introduction 

Ultra-relativistic collisions of nuclei, so called “Little Bangs” are observed at the Relativistic 
Heavy Ion Collider (RHIC) of the Brookhaven National Laboratory, New York. The PHENIX 
Experiment at RHIC detects and identifies a large variety of the outcoming particles over a broad 
energy range. PHENIX has been participating in the Beam Energy Scan program of RHIC, and 
collision data of many nuclei (d. He, Cu, Au, U) were taken at several center of mass energies (per 
nucleon), from 7.7 GeV to 200 GeV. The goal behind this variety of energies is to investigate the 
phase structure of QCD matter and the nature of the quark-hadron transition. This transition is a 
cross-over at low baryochemical potential [1], corresponding to top RHIC energy. It is however 
believed to be a first order phase transition at higher baryochemical potential (i.e. lower collision 
energies), while there may be a critical endpoint (CEP) somewhere in between. Correlation and 
fluctuation measurements [2, 3] are a powerful tool to be utilized in the search for direct evidence of 
the phase transition and the location of the critical point. In the present paper we focus on the soft 
pion production in these collisions, in particular the measurement of Bose-Einstein correlations 
between these low momentum pions. We discuss general features of Bose-Einstein correlations 
below, and then show results of various Bose-Einstein correlation measurements. 

2. Basics of Bose-Einstein correlations and HBT radii 

R. H. Brown and R. Q. Twiss, radio astronomers, used two radio antennas pointed at a distant 
star and chose to detect each signal separately to take away its high frequency component. They 
saw an interference term present in the product of these two detected signals, and it turned out, 
that this interference enabled them to measure the angular radius of the observed star [4]. This 
technique was then translated to high energy physics by Goldhaber et al. [5], and it turned out, that 
Bose-Einstein correlations of pion pairs coming out of high energy nuclear collisions can be used 
to study the space- time extent and shape of the particle emitting source. 

The basic idea can be described in a simple quantum-mechanical picture: the detected two- 
particle momentum distribution depends on the two-particle probability amplitude, i.e. the absolute 
square of the two-particle wave function. This wave function has to be symmetric for the exchange 
of the two particles, and in the end this results in a simple connection between the normalized 
particle emitting source S{r) and the momentum-correlation function C{q)'. 

C(^) = l + |%)|' (2.1) 

where S{q) is the Eourier transform of S{r), and q = p\ — p 2 the four- momentum difference 
of the two identical bosons. In case of usual sources, the C(^) —)• I for q ^ oo, and C(0) =2 (as 
5(0) = I for normalized sources). 

There are however many complications that distort the above simple picture. One of these is 
the Coulomb-interaction (and other final-state interactions, which however can be neglected in case 
of pions) between the two particles. This can be handled by using the solution of the two-paiticle 
Schrodinger-equation with a Coulomb-potential, as discussed in Ref. [6] for details. Another com¬ 
plication is, that a significant fraction of the produced pions comes from the decays of long lived 
resonances. These pions are produced very far from the hot, hydrodynamically expanding core, and 
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give a very narrow contribution to C{q). This narrow component cannot be resolved experimen¬ 
tally, thus the correlation function seems to have an intercept smaller than 2. Such a two-component 
source can be described in the core-halo picture (introduced in Refs. [7, 8]), where the experimen¬ 
tally resolvable correlation function is 

C(^) = l+A|Seore(^)|^ (2-2) 

where Score is the normalized source function of primordial pions, 0 < A < 1, and C(g —)• 0) = 1 -|- A 
for the experimentally resolvable q region. 

Usually, the ID correlation function is measured as a function of ^inv = y/—q^, but 3D infor¬ 
mation on the source can also be extracted from measurements, if statistics permits it. The z or 
longitudinal {long) direction is that of the beam. The plane perpendicular to it is parametrized with 
a pair-coordinate system introduced in Refs. [9, 10]. In this system the out direction is that of the 
average transverse momentum of the pair, while the side direction is perpendicular to both out and 
long. Also, instead of the center of mass or the lab system, the longitudinally comoving system 
(LCMS) is used usually (introduced in Ref. [11]). In a Bertsch-Pratt LCMS system, the average 
momentum is while from the mass-shell condition one obtains 


qo — ,. 9out — Ptqout 

Mf 


(2.3) 


i.e. the time component of the momentum difference is related to the difference in the outward 
direction. 

The main physics information is to be extracted from the shape of the measured correlation 
functions in one dimension (as a function of ^inv) or in three dimensions (as a function of ^out> 9side 
and ^long)- The simplest case is to measure the width’ of the correlation function (the so-called 
HBT radii named after H. Brown and Twiss), with a Gaussian assumption. This is not always valid 
however, and in an ideal case, model independent shape analysis has to be pursued [12]. Others 
restore the source image from the correlation function, via the method of imaging [13]. In this 
paper we report about results of ordinary Gaussian analyses, and the measurement of Gaussian 
HBT radii. 

And how are HBT radii connected to the the critical point or phase transitions? It turns out, 
that in a simple hydrodynamic picture, the out and side directed HBT radii are [14] 


p2 _ 

^out 


l + ^uj 




and 


p2 _ 

^side 


i+t«r 


(2.4) 


which implies two things: radii scale as 1 /y/in^ with the average transverse mass of the pair, and 
out minus side difference is related to the emission duration At. This is particularly interesting, 
since a strong first order phase transition would cause a long emission and thus a large difference 
of these radii - which was not the case experimentally [15]. Readers shall however be reminded, 
that there is rich physics hidden in this story, and are encouraged to read Ref. [16] for an in-depth 
analysis. 


3. PHENIX HBT results 

In this section, we present recent PHENIX results on Bose-Einstein correlations. We show a 
few example correlation functions in Eig. 1, from Ref. [17]. These (and others) were then analyzed 
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Figure 1: A variety of three dimensional correlation functions at a glance, from Ref [17]. 


by extracting their Gaussian width, i.e. the HBT radii. In the following subsections, we present 
such results on HBT radii. 

3.1 Transverse momentum scaling 

Correlation functions of pion pairs of several average transverse mass bins were measured in 
39 to 200 GeV Au+Au collisions [17]. The nit dependence of the resulting Gaussian radii /?out^ 
/?side> and 7?long^ compared to recent STAR data [18] are shown in Fig. 2. The data of different 
centrality selections show characteristic 1 /^JWt scaling patterns as a function of the transverse mass 
nit of the emitted pion pairs, consistent with hydrodynamic expansion. A good level of agreement 
is seen between the PHENIX and STAR data sets, while the PHENIX data extends the covered ni, 
range of RHIC HBT radii. It is important to note that the scaling also holds in 200 GeV d+Au 
collisions [23]. 

Correlation functions of kaon pairs were also recently measured at PHENIX. Since the pion 
sample is more affected by hadronic rescattering and the decays of long lived resonances, kaons 
represent a more clean proble of the particle emitting source. PHENIX recently measured kaon 
HBT radii at several centralities of 200 GeV Au+Au collisions [19]. Eig. 3 shows a comparison 
of the nit dependence of HBT radii between charged pions and kaons. There is no significant 
difference in the side direction (i.e. scaling of different particle species holds), while the kaon 
radii in the out and longitudinal directions are slightly larger than the pions (at the same nit). The 
difference is more pronounced in central collisions. These results may point to different dynamics 
and duration of kaon production as compared to pion production. 

3.2 Azimuthally dependent HBT 

HBT radii were recently measured at PHENIX with respect to the 2nd- and 3rd-order event 
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Figure 2: Transverse momentum scaling of HBT radii from 39 to 200 GeV, from Ref. [17]. 
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Figure 3: Transverse momentum scaling of kaon radii, from Ref. [19]. Significant difference from pion 
radii is observed in central collisions, in the out and longitudinal directions. 
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Figure 4: Azimuthally sensitive HBT radii w.r.t. the 2nd and 3rd order event planes (top panel), from 
Ref. [20]. Their first Fourier-components as a function of initial eccentricity are shown in the bottom panel. 


plane in 200 GeV Au+Au collisions [20]. A strong azimuthal-angle dependence of the Gaussian 
HBT radii was observed with respect to both the 2nd- and 3rd-order event planes, as shown in Fig. 4 
(toppanel). The oscillation amplitudes were extracted by fitting 7?^ q -|- 2 ^ 7?^ „cos[?i(0 — 'F„j)] 

n—m^2m 

type of functional forms. The amplitude ratios versus initial spatial anisotropies £„ 

(from Monte-Carlo Glauber simulations) are shown also in Fig. 4 (bottom panel). Since these 
ratios are sensitive to the final eccentricity, the results for n = 2 indicate that the initial eccentricity 
is reduced during the medium evolution, but is not reversed. Results for n = 3 however indicate, 
that the initial triangular shape may be reversed by the end of the medium evolution, and that 
the 3rd-order oscillations are largely dominated by the dynamical effects from triangular flow, as 
detailed in Ref. [20] 

A new technique in measuring azimutally sensitive HBT radii is the application of event shape 
selection. In this method one controls not just the direction of the event plane in each event, but 
also the magnitude of the second order flow vector Q 2 , measured by the Reaction Plane Detector 
(RxNP). Q 2 is defined as ^;W;(cos(2(/),),sin(20,)) with w, being the multiplicity in the /-th segment 
of RxNP. The measured eccentricity (defined as 27?^ „/7?y q) was found to be enhanced for events 
with higher magnitude of the flow vector Q 2 , as detailed in e.g. Ref. [21]. 

3.3 Geometric scaling 

1 /3 

As the initial system volume scales with Apart, the dependence of the HBT radii is often 
analyzed. Recently, HBT radii measurements in 39, 62 and 200 GeV Au-i-Au collisions were com¬ 
pared [22], as shown in Fig. 5 (top panel). All energies show a linear dependence, approximately 
independently of the beam energy. This indicates that at these RHIC energies, the dynamics of the 
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Figure 5: System size scaling from Refs. [22, 23]. Linear scaling is observed, but d+Au data seem to fall 
on a different scaling curve. 


system does not change dramatically, similar Apart means similar initial system size, which in turn 
leads to a similar homogeneity length. 

PHENIX also measured three dimensional HBT radii in d+Au collisions as a function of cen- 
trality [23]. The linear Ap'rt dependence shows a similar slope in Au+Au and d+Au collisions (at 
the same average transverse momentum). The similarity of the slope is stronger in case of /?out and 
f^side (hinting at similar transverse dynamics of the two systems), while the /?iong slope is slightly 
different, maybe due to different longitudinal dynamics in d+Au and Au+Au collisions. Results 
are shown in Fig. 5 (bottom panel). 

Besides the above system size scaling, more detailed geometrical scalings were discovered. 
These are based on the characteristic initial transverse size, defined as 1//?^ = 1/a^ + l/a^ with 
Ox,y being the RMS widths of participant nucleon distributions in the x and y directions (estimated 
via Glauber Monte Carlo calculations). The space-time extent of the emission source at freeze-out, 
as measured by the HBT radii, reflects the initial size of the system, the expansion time or lifetime 
T, as well as transverse expansion rate (visible via the \/x/tnt dependence of the HBT radii, as 
discussed in Section. 2). Since the expansion time is expected to scale with R, the HBT radii are 
also expected to scale with R in a given nit range. This scaling was indeed observed at PHENIX [17] 
as shown in Fig. 6. Note furthermore, that it also holds in case of d+Au collisions [23]. The detailed 
scaling laws indicate that R is a better predictor of final state homogeneity length, than Apa^. which 
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Figure 6: Geometric scaling versus R, from Ref. [17]. 
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Figure 7: Excitation function of different combinations of HBT radii, from Ref. [17]. Non-monotonicity is 
seen around .Vsnn =40 GeV. 


is consistent with observations that HBT is sensitive to the expansion dynamics. 


3.4 Excitation function 

In order to access physical descriptors of the expanding medium like emission duration and 
expansion velocity, quantities (^side — V2^)/^iong were measured, as discussed 

in Ref. [17]. These were investigated as a function of collision energy, and results at nit = 0.26 
GeV/c^, for 0-5% centrality are shown in Fig. 7. The observed non-monotonic pattern near 
^5 nn =40 GeV could be an indication of the vicinity of the softest point of the QCD EoS. It 
is however important to note that the measured HBT radii are smaller than the actual sizes of the 
system, and are linked to expansion dynamics. Due to the amount and depth of physics involved in 
these measurements, further experimental studies and full hybrid models are required to confirm if 
the observed non-monotonic patterns are linked to the critical end point in the QCD phase diagram. 
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3.5 Levy exponents 

At the critical point universal power law behavior of several observables can be seen as a 
function of the temperature of the system, characterized by critical exponents, one being that of 
the spatial correlation function. Levy distributions are present in the experimental momentum 
correlation functions and the spatial correlation functions restored from it. The Levy-exponent 
a of these correlation functions is equivalent to the spatial correlation exponent, conjectured to 
be near 0.5 in case of the critical end point [24]. PHENIX performed measurement of the Levy- 
exponent in 200 GeV Au-i-Au collisions [25]. At this energy, the exponent is far from both the 
Gaussian case (a = 2) and the critical case (a = 0.5). A detailed measurement of the excitation 
function of the correlation exponent is expected to substantially refine results on the QCD critical 
end point. 

4. Beyond HBT 

As present conference volume is about the search for and characterization of the QCD critical 
end point, this paper cannot be complete without mentioning some more related PHENIX results, 
beyond HBT. 

One set of important results is that of neutral pion spectra, which were measured in 39 GeV to 
200 GeV Au-i-Au collisions. Midrapidity nuclear modification factors measured in the 1 < p, < 10 
GeV/c range were compared to each other in Ref. [26]. Eor central collisions, suppression is 
significant at y^^NN =39, 62 and 200 GeV as well. On the other hand, in mid-peripheral 39 GeV 
collisions, the nuclear modification factor Raa is consistent with unity above pt >'i GeV/c. 

PHENIX recently measured the anisotropic flow coefficients v„ for identified pions, kaons and 
profons, relative to the n-th order event planes, in 200 GeV Au-i-Au collisions [27]. These flow 
coefficients show characteristic patterns consistent with hydrodynamical expansion of the matter 
produced in the collisions. They all fall on modified quark number scaling curves, if v„/{nq)"/'^ is 
plotted versus KE^ jnq. This may indicate a hydrodynamic expansion, a quark coalescence process, 
and/or an acoustic nature of the QGP Interestingly, these scalings with the number of constituent 
quarks hold for n = 2,3 also in 39 and 62 GeV Au-i-Au collisions [28]. 

5. Summary 

A wide variety of HBT data is now available. Geometric and hydrodynamic scalings are 
observed for a broad y^^NN range: HBT radii versus 1 / versus Vp^-t and versus R follow 
simple linear curves. Excitation function of HBT radii however show a non-monotonic behavior 
near = 40 GeV. This may be attributed to non-monotonicity in the emission duration; more 
detailed analysis and measurements have however to be performed to understand and confirm the 
relation of this behavior to the QCD critical end point. Measurements of the correlation exponent 
could for example substantially refine our understanding. There are however other probes measured 
in the vicinity of this point. Non-monotonic behavior is seen in neutral pion nuclear modifications, 
while quark number scaling is preserved from 200 GeV to 39 GeV. All of these observables and 
their excitation functions shuld be understood in a global complex picture, in order to acquire 
detailed knowledge of the QCD critical end point. 
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